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We show that self-localized ground states can be created in the spin-balanced gas of fermions 
with repulsion between the spin components, whose strength grows from the center to periphery, in 
combination with the harmonic-oscillator (HO) trapping potential acting in one or two transverse 
directions. We also consider the ground state in the non-interacting Fermi gas under the action of 
the spatially growing tightness of the one- or two-dimensional (ID or 2D) HO confinement. These 
settings are considered in the framework of the Thomas-Fermi-von Weizsacker (TF-vW) density 
functional. It is found that the vW correction to the simple TF approximation (the gradient term) 
is nearly negligible in all situations. The properties of the ground state under the action of the 2D 
and ID HO confinement with the tightness growing in the transverse directions is investigated too 
for the Bose-Einstein condensate (BEC) with the self-repulsive nonlinearity. 
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I. INTRODUCTION 

One of fundamental conclusions produced by ongoing 
studies of the dynamics of matter waves in Bose-Einstein 
condensates (BECs) is that the effective local nonlinear- 
ity, induced by inter- atomic collisions, gives rise to robust 
bright solitons and soliton complexes [lH7|. In addition 
to their fundamental significance, the solitons may be em- 
ployed in applications. In particular, the use of solitons in 
matter-wave interferometers should help to dramatically 
increase the accurac y o f these devices, see Refs. (8|-[Tl| 
and a recent review [13|. 

In the usual settings, the existence of bright solitons re- 
quires the presence of the self-focusing nonlinearity. They 
may also be supported by nonlinear lattices, which fea- 
ture alternation of spatial domains of self-focusing and 
defocusing, which give rise to the corresponding peri- 
odic pseudo-potential 0, [ill- On the other hand, self- 
defocusing nonlinearities, acting in the combination with 
periodic linear (lattice) potentials, support bright soli- 
tons of the bandgap type [1, H, Nevertheless, a com- 
mon belief was that the self-defocusing per se could not 
give rise to bright solitons. The situation had changed 
when it was demonstrated that the repulsive cubic non- 
linearity with the local strength growing with the dis- 
tance from the center, r, faster than r^, where D is the 
spatial dimension, can readily support stable solitons and 
solitary vortices (TBI. [l6j. In the same works, it was pro- 
posed how the corresponding profiles of the nonlinearity 
modulation can be created in BEC and nonlinear optics. 
In particular, the spatial modulation of the scattering 
length of inter-atomic collisions, induced by spatially in- 
homogeneous magnetic or optical fields via the Feshbach- 
resonance mechanism, can give rise to the required profile 
in BEC. Similarly, robust solitons and vortices were pre- 
dicted in a more exotic model, with the constant strength 



of the defocusing nonlinearity but the diffraction coeffi- 
cient decaying faster than WA - 

It may be interesting to implement this new option for 
the creation of bright solitons by purely repulsive nonlin- 
earity landscapes in other physical settings. For instance, 
it was recently demonstrated that the same mechanism 
works too in one- and two-dimensional (ID and 2D) me- 
dia with the repulsive quintic nonlinearity [l8| . 

The present work aims to propose a way of creating 
bright solitons in Fermi gases, with balanced spin-up and 
spin-down components. The possibility is to combine the 
spatial growth of the s-wave scattering length, which ac- 
counts for the repulsion between the components, in one 
or two directions, and the harmonic-oscillator (HO) con- 
fining potential acting in the other directions. It should 
be noted that the dilute two-spin-component fermionic 
gases with the repulsive s-wave interaction do not feature 
phase coherence, as they are not superfiuids [l^. Never- 
theless, their ground states can be accurately described 
by means of the density- functional theory j20l - |22| . 

Here we adopt the Thomas-Fermi-von Weizsacker (TF- 
vW) form of the density functional to predict density pro- 
files and chemical potentials of the Fermi droplets (i.e., 
bright solitons) in two different configurations, which 
are considered in Sections II and HI, respectively: the 
quadratic growth of the scattering length in one direction 
(z) and transverse HO confinement in the (x, y) plane; 
or the quadratic growth of the scattering length in the 
{x,y) plane and HO confinement acting along z. Then, 
in Sections IV and V, following previous results obtained 
in the framework of the mean-field description of BEC 
j23l - l25j . we consider Fermi gases without intrinsic inter- 
actions (in particular, it may be a spin-polarized, i.e., 
single-spin-component, gas) and demonstrate that they 
give rise to localized ground states under the action of the 
2D or ID HO confinement, if its tightness, i.e., the corre- 
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spending confinement frequency, grows in the transverse 



directions faster than 



respectively. In addi- 



tion, we demonstrate that the same confinements with 
the spatiahy growing tightness support, in a similar way, 
localized ground states in BEC with the self-repulsive 
nonlinearity. 



II. ID NONLINEARITY MODULATION WITH 
THE 2D TRANSVERSE CONFINEMENT FOR 
THE SPIN-BALANCED INTERACTING FERMI 
GAS 

In the presence of an external trappin g p otential 
C4xt(a;, y, z), the Hohenberg-Kohn theorem [201 ensures 
that the single-body local density y, z) of the ground 
state of a quantum system composed of interacting iden- 
tical particles can be obtained by minimizing the energy 
functional, 



E\n] 



F[n\ + J J J Uc^tix,y,z) n{x,y,z) dxdydz, 

where F[n] is an internal-energy functional, which is in- 
dependent of Ucyit{x,y, z). For the dilute normal Fermi 
gas with two equally-populated spin states, the TF-vW 
internal energy at zero-temperature is written as [l9l | 



F\n\ 



5 2m 8m n 

dx dy dz, 



(2) 



with the nonlinearity strength determined by the s-wave 
fermion-up-fermion-down scattering length, which (as we 
assume in this work) may be modulated along the axial 
direction, a^^{z): 



g{z) = [ATrh^/m) afi{z). 



(3) 



The first term in functional is the TF kinetic energy 
of the Fermi gas at zero temperature, while the second 
term is the vW correction to the kinetic energy of the 
inhomogcneous gas (the surface term). For the coefficient 
in front of the vW term, the phenomenological value, A = 
1/3, may be adopted for non-superfluid fermions j2l| . 
The third term in ^ is the mean-field approximation 
of the s-wave interaction between fermions with opposite 
spins. 

The HO confinement in the transverse plane is ac- 
counted for by the external potential, 

U,^tix,y) ^ {l/2)mL,lr\ (4) 

with r — {x'^ + y^Y^"^ and confinement radius aj_ = 
\Jh,l (mwj,). We suppose that the spatially modulated 
scattering length also features the quadratic coordinate 
dependence, with a characteristic scale ao > 0, 



(5) 



By minimizing the energy functional ([2]), which is sub- 
ject to the normalization constraint, imposed by the fixed 
total number of fermions. 



N = 



7i(x, y, z) dxdydz, 



(6) 



one derives the equation for the local density, 

-AV^ + (37r2)2/3n2/3 ^ ^2 ^ ^^2 1 ^ ^ ^^y^, (7) 

with chemical potential /i fixed by normalization 
Equation ([7]) is written in the scaled notation, with 
lengths measured in units of aj_ and energies in units 
of huj±. In this form, the equations depend on the single 
parameter, viz., the adimensional interaction strength. 



7 = 47rao/a^. 



(8) 



Multiplying the left-hand side of Eq.© by y/n and inte- 
grating the result over the spatial coordinates, we obtain 
an expression for the chemical potential in terms of den- 
sity n{x, y, z): 



1 
2N 



+r'^n + 7Z^n'^] dx dy dz. 



-AV^V^ (V^) + (37r2) 



2^2/3 ^5/3 (9) 



(10) 



In the TF regime, when the vW correction is negligible, 
Eq. ([7]) for the density profile reduces to a cubic algebraic 
equation for v}/^, 

(37r2)2/^i2/3 + +7^2^ = 2/i, (11) 

at r > ttf, where the TF radius 



at r < rxF, and n 
is 



2^. 



(12) 



To find solutions of full equation ([7]), including the 
vW correction, an imaginary-time-derivative term was 
added to it, and ensuing stationary solutions were ob- 
tained by means of the finite-difference Crank-Nicolson 
predictor-corrector method, which keeps the fixed value 
of N [2^ . In Fig. [1] we plot the so obtained 3D den- 
sity profile n(r, z) of the spin-balanced Fermi gas, pro- 
duced by the numerical solution of Eq. ([7]) for /Lt = 10, 
the nonlinearity strength 7 = 1, and A = 1/3 (as said 
above) , which corresponds to the self-trapped mode built 
of = 612 fermions. The figure shows that, while 
along the r direction the density practically vanishes at 
ttf = \/20 ~ 4.47, as predicted by the TF approxima- 
tion, along the z direction the density vanishes only at 

\z\ — !• 00. 

In Fig. [2] we plot the radial density, n{r,z = 0), for 
N = 100 fermions, and compare the solutions of the 
full equation ([7]) with the TF approximation based on 
Ea. ([TT|) . The figure shows that the vW gradient term 
mainly affects the behavior of the density near the sur- 
face layer (see the insets in the figure): instead of vanish- 
ing at a finite distance rxF from the center, the density 
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FIG. 1: (Color online). The three-dimensional density profile 
n(r, z) of the two-component spin-balanced Fermi gas, under 
the action of the ID nonlinearity modulation, see Eq. ([5]), and 
the harmonic-oscillator confinement in the transverse plane 
[Eq. Q]. The chemical potential is /i = 10, adimensional 
nonlinearity strength j — 1 [see Eq. ((S))], and A = 1/3. 



FIG. 3: (Color online). Chemical potential /i (in units of fiij_L) 
vs. the number of particles, A'^, for the spin-balanced Fermi 
gas with nonlinearity strengths 7 = 5, 1 and 0.5, in the same 
settings as in Figs. [l]and[2l Lines represent the dependences 
obtained from the full Thomas-Fermi-von Weizsacker model, 
i.e., Eq. ((Tjl with A = 1/3, while points correspond to the 
TF approximation (A — 0). The inset is a blowup of the 
dependence at small values of A*". 




FIG. 2: (Color online). Radial density n(r, 0) of the spin- 
balanced Fermi gas with nonlinearity strength 7 = 1 and 
number of particles A'' = 100, in the same setting as in[T] Solid 
lines: solutions of the full equation (O, which includes the 
von Weizsacker term with A — 1/3. Dashed lines: solutions 
produced by the Thomas-Fermi (TF) approximation based on 
Eq. (fTT|) [i.e., Eq. ^ with A = 0]. Insets display a blowup of 
the structure of each state across the surface layer. 



repulsive nonlinearity. This expectation is confirmed by 
Fig. 21 In the upper and lower panels of this figure, we 
plot the radial and axial density profiles, viz., n{r, z = 0) 
and n{r — 0,z) respectively, for N = 100, A = 1/3 (i.e., 
the vW is taken into account here, although it produces 
very little change) and three values of the nonlinearity 
strength, 7 = 0.5, 2, and 5. 

As said above, chemical potential ^ determines the to- 
tal number N of atoms in the ground state, hence its 
width increases with /i, for fixed nonlinearity strength 7. 
In Fig. Owe plot the root-mean-square (rms) radial size 
of the ground state, (r^)^/^, versus N, for three values of 
the nonlinearity strengths, 7 = 5, 1 and 0.5, using the 
the full equation ([7]) with A = 1/3 (lines in Fig. [S]) and 
the TF approximation corresponding to A = (points in 
Fig. [5]). The results produced by the full and simplified 
models are essentially the same for different values of 7 
and, virtually, for all N. 

III. 2D NONLINEARITY MODULATION WITH 
THE ID TRANSVERSE CONFINEMENT 



vanishes at r — > 00. Nevertheless, this effect is weak, and 
it becomes negligible for larger N. 

Similar features are exhibited by the dependence of 
the chemical potential, fi, on the number of atoms, N, as 
shown in Fig. [S] The change of these dependences caused 
by the vW term is very small even for the state built of 
a dozen of atoms, see the inset in the figure. Thus, the 
simplified TF functional, disregarding the vW corrections 
(A = 0), is sufficient for the study of the ground state in 
the present model. 

Regarding the ground state, it is natural to expect that 
its width reduces with the growth of strength 7 of the self- 



We now consider the setting with the HO confinement 
acting along the longitudinal direction, viz., 

= il/2)mujy, (13) 

cf. Eq. (j4]). On the other hand, the spatial modulation 
of the scattering length, a^i, is adopted here to be two- 
dimensional, cf. Eq. ([S]): 

auir) = (ao/al) r\ (14) 

with flo > 0, where = yjh/ (niLu^) is the characteristic 
length of the longitudinal HO confinement (|13p . Thus, 
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FIG. 4: (Color online), (a) Radial n{r,z — 0) and (b) axial 
7i(r — 0, z) density profiles of the ground state composed of 
A'^ — 100 fermions, for A = 1/3 in Eq. ([7]) and three values of 
nonlinearity strength 7, which is defined as per Eq. ([8}. The 
setting is the same as in Figs. [l]|3l 
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FIG. 5: (Color online). The rms radial size {r^}^^^ vs. the 
number of fermions, A*", for nonlinearity strengths 7 = 5, 1 
and 0.5. Lines are produced by full equation (O with A = 1/3, 
while points correspond to the TF approximation, A = 0. 
The inset shows the magnification of the dependence at small 
values of A''. The setting is the same as in Figs. IH4I 



the present setting is a reverse of that considered in the 
previous section. 



The corresponding TF-vW internal energy is 

"3 ft2 ^2 



F[n] 



^ (37r2)2/3n5/3 + A^^^15) 
5 2to^ ^ 8m n ^ ' 



dx dy dz, 



(16) 



with g{r) given by Eq. ([3]), in which a^\,{z) is replaced by 
a-|-j,(r), see Eq. p4)) . By minimizing the fuh energy func- 
tional, E\n] = F[n] + J J J Ucxt{z)n{x,y, z)dxdydz, with 
the constraint of the normalization of n(r), we obtain 



-AV^ + (37r2)2/3n2/3 + + Tr^n] = 2fiV^, (17) 



cf. Eq. d?]). In Eq. ([T7)l lengths are measured in units 
of az, and energies in units of hujz, which leaves, as the 
single control parameter, the adimensional strength of 
the nonlinearity, 



r = Anai/az 



(18) 



cf. Eq. (El). 

In the TF approximation, which, as before, neglects 
the vW gradient term, the density profile, n{r,z), satis- 
fies an algebraic equation [cf. Eq. (fTT|) ]: 



(37r2)2/3n2/3 + z^ + Tr^n = 2^. 



(19) 



Setting n{r, z) = in Ea.([T^. one gets the TF axial size 
of the ground state, ztf = \/^^ with 71 = at \z\ > ztf 
in the TF approximation, cf. Eq. (fT2l). 

In Fig. [6] the 3D density profile n{r, z) of the Fermi 
gas is plotted for /i = 10, F = 2, and A = 1/3 in Eq. 
(fT7| . The figure shows that, while along the z direction 
the density practically vanishes at ztf = V20 ~ 4.47, 
along r the density vanishes only at r — > 00. 

In Fig. [7]we plot axial density n(0, z) for two values of 
the chemical potential, fi ~ 10 and 20, in panels (a) and 
(b), respectively, both taken with F = 2. In this figure 
we compare solutions of the full equation p7)) to their 
counterparts produced by the TF approximation, which 
is based on Eq. ([T^ . The figure shows, once again, that 
the vW gradient term mainly affects the behavior of the 
density in the surface layer: Instead of vanishing at a 
finite distance (ztf) from the center of the cloud, it van- 
ishes at \z\ — 00. However, this effect is weak, becoming 
negligible at larger values of the chemical potential. 

As shown in Fig. |51 a similar phenomenon is observed 
in the dependences of the chemical potential fi on N, 
which were produced by the full equation [T7] with A = 
1/3 (lines in Fig. |5]), and by the TF approximation 
(points in Fig. |5]). The effect of the vW gradient term 
remains inconspicuous even at small values of N, and at 
all values of F. 

A natural expectation that the radial width of the 
ground state reduces with the increase of the nonlinear- 
ity strength, F [i.e., the gas is stronger compressed by 
the effective nonlinear (pseudo)potential], is confirmed 
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FIG. 6: (Color online) Three-dimensional density profile 
n(r, z) of the two-component spin-balanced Fermi gas under 
the ID external confinement, see Eq. (|13p . and 2D modu- 
lation of the local nonlinearity in the transverse plane [Eq. 
(|14|l ]. The chemical potential is /i = 10, adimensional nonlin- 
earity strength is F = 2, and A = 1/3 is adopted in Eq. (|17|l . 
Lengths are measured in units of a^, and energies in units of 
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FIG. 8: (Color online). Chemical potential jj, versus number 
of fermionic atoms for the two-component spin-balanced 
Fermi gas under the ID confinement and 2D modulation of 
the local nonlinearity. The nonlinearity strength takes values 
F = 5, 2 and 0.5. Lines are obtained from full equation (|17|l . 
while points correspond to the TF approximation (A — 0). 
Energies are in units of fkuz. 




FIG. 7: (Color online). Axial density profiles n{0,z) of the 
two-component spin-balanced Fermi gas under the ID con- 
finement [Eq. (|13|) ] and 2D modulation of the nonlinearity 
[Eq. (O], whose global strength is F = 2, [see Eq. (fT5)) . for 
two values of the chemical potential: fi — 10 (a) and = 20 
(b). Solid lines: solutions of Eq. (|17|l with A = 1/3; dashed 
lines: solutions of Eq. (fT9)) . i.e., Eq. (fTfl) with A = 0. Units 
are the same as in Fig. [G] 
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FIG. 9: (Color online). The rms size (z^)^^'^ vs. number A'^ of 
fermions of the ground state with nonlinearity strengths F = 
5, 2 and 0.5. Lines and points represent the results obtained, 
respectively, from the full equation (|17p with A — 1/3, and 
from the TF approximation corresponding to A = 0. The 
inset shows details of the dependences at smaller values of A''. 
The setting is the same as in Figs. [6][8l 



by Fig. |9l In this figure we plot the rms size of the 
ground state, (z^)^/^, versus the number of particles, N, 
for three values of the nonlinearity strength, F = 5, 2 and 
0.5, using the fuU equation ^ with A = 1/3, the TF 
approximation corresponding to A = (lines and points, 
respectively, in Fig. Thus, we conclude that the effect 
of the vW term is inessential in this setting too. 
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IV. SELF-TRAPPING IN THE 
NON-INTERACTING FERMI GAS AND 
SELF-REPULSIVE BEG DUE TO THE AXIAL 
MODULATION OF THE 2D CONFINEMENT 

A. The non-interacting Fermi gas under the 2D 
confinement with the strength growing in the axial 
direction 

In the spin-balanced ideal Fermi gas, without any in- 
teraction between the two spin components, or the single- 
component spin-polarized gas, in which direct interac- 
tions are suppressed by the Pauli principle, it is possi- 
ble to induce self-trapping along the longitudinal direc- 
tion (z) by introducing a z-dependent modulation of the 
transverse HO confinement frequency, u}±. The corre- 
sponding energy functional is 



E[n] = 



3 ft2 
5 2m' 



ft2 (Vn)' 



(20) 



dx dy dz 



(21) 



By minimizing this functional, subject, as before, to the 
constraint of the normalization of n{r), one obtains 



(Stt' 



)2/3^2/3 



2fiy/n, (22) 



where ^ is the chemical potential fixed by the total num- 
ber of atoms, N [see Eq. ©]. Solving Eg. d^^ . we aim 
to produce the longitudinal density profile, 



ni{z) 



n(r) dxdy. 



(23) 



The TF approximation corresponds, as above, to set- 
ting A = in Eq. ([22| , which makes it a simple algebraic 
equation that immediately yields an explicit solution [this 
was not available in the presence of the interaction be- 
tween the spin components, cf. Eqs. ([7|) and (fTTj)]: 



n (r, z) 



^ f (3^2-)-! _ tolizyf'^ , at Lol{zy < 2fi, 
\ 0, at a;^(z)r2 > 2^ 

(24) 

(obviously, the solution exists only for > 0). This 
solution is physically meaningful if its norm converges, 
which implies the convergence of r^p (z') dz' = 
2fj, ujj^ {z') dz' at z — >■ oo. The eventual condition 
is that ci;^(z) must grow faster than |z|, i.e.. 



uj±{z)/ 



z — )• CX), 



(25) 



at |z| oo. Further, the substitution of explicit solution 
([M)) into Eq. makes it possible to obtain the TF 

approximation for ni in an explicit form too: 



ni(z) 



16^2 (mii\5n 4 



(26) 




FIG. 10: (Color online). Profiles of the longitudinal density 
m {z) of the non- interacting Fermi gas, with the z-dependent 
modulation of the transverse-confinement frequency, chosen 
as per Eq. (|27|) with a = 0.01 and 13 = 1, for three values of 
the chemical potential: /i = 8, 6.5, and 5. Lines: a numerical 
solution of Eq. (|22p for A = 0. Points: analytical result 



Below we consider the most natural modulation form, 

a;i(z) ==az2+/3, (27) 

with positive a and (3. Obviously, it satisfies condition 
In Fig. [TU]wc plot longitudinal profiles of the re- 
duced density (|23p for the z-dcpendent modulation of the 
transverse-confinement strength given by Eq. (|27p with 
a = 0.01 and (3 = 1. Dots in Fig. [TU] represent analytical 
result (pS)) , and lines depict (for the sake of checking the 
correctness of the analytical result) the numerical solu- 
tion of Eq. (gl]) with A = 0, obtained as in Ref. ^ [in 
the latter case, the 3D density is numerically integrated 
in the transverse plane to reduce it to ni, see Eq. ([^5]). 

To address effects of the vW term in the present case, 
in Fig. [TT] we plot chemical potential /i as a function of 
the number of atoms N, using the numerical solutions of 
Eq. dm with A = 1/3 and A = 0. It is seen that the 
effect of the vW gradient term is again very weak, for any 
number of particles. For the same purpose, in Fig. I12l we 
plot the radial density, n(r, z = 0), as found from the full 
equation and produced by the TF approximation, for the 
same set of values of the chemical potential as in Fig. (TU] 



B. A self-repulsive bosonic condensate under the 
2D confinement with the strength growing in the 
axial direction 

A mechanism of the formation of the localized ground 
state, similar to that presented in the previous section, 
may be applied to BEG with repulsive interactions be- 
tween atoms, subject to the action of the 2D confinement 
in the plane of (x, y), with the tightness growing along z. 
The respective effective ID equation for the mean-field 
wave function, $ (z,t), was derived in Ref. [1^, starting 
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FIG. 11: (Color online). Chemical potential /i vs. the num- 
ber of fermionic atoms, A^, of the ground state of the non- 
interacting Fermi gas subject to modulation (|27p of the fre- 
quency of the transverse confinement. As indicated in the 
figures, the TF approximation (A = 0) produces the depen- 
dence which is very close to that generated by the full model, 
including the von Weizsacker term (A — 1/3). The inset is a 
blowup of the dependence at small values of A'^. 



(in Ref. [25|, the norm of the ID wave function was 1, 
while as in Eq. ([28| was replaced by Nag). Looking for 

— 1/2 

a stationary state as $(z,i) = Os cxp{~ifj,T)U{Z), 
<t and Z = z/as, we arrive at an 



where T = {h/r 
equation for real function U{Z): 



4C/2f7, 



(30) 



with 57 = (mal/fij uj±. The TF approximation for so- 
lutions of Eq. ([50]) is obvious: 



. . ^ / {2n^iz)) ^^i'^-niiz), at n^{z) < /i, 
^ ' I 0, at n^{z) > n, 

(31) 

cf. Eq. dMl)- 

If, in particular, io±{z) grows faster than \/\z\ at \z\ — 
oo [i.e., it obeys condition (|25p ]. the TF approximation 
gives rise to a simple asymptotic dependence of the total 
number of bosonic atoms on the chemical potential: as 
follows from the substitution of expression (j3ip into Eq. 



A,=1/3: n=8 
1^=6.5 
|j,=3 
X=0: [1=8 
|a=6.5 
H=3 




p=1 



FIG. 12: (Color online). Radial density n{r,z — 0) of the 
non-interacting Fermi gas, with the modulation of the fre- 
quency of the transverse confinement defined by Eq. (|27p . 
for three values of the chemical potential: ^ — 8, 6.5 and 3. 
As above, symbols A = 1/3 and A = designate, severally, 
the results produced by the full equation (|22|l , and by the TF 
approximaiton (|24p . Units the same as in Fig. 1101 



/+00 2 r+oo 

U'iZ)dZ^i^ . 



dZ 



nliz) 



(32) 



for n large enough. On the other hand, if ujj_ grows slower 
at \z\ oo, namely, n^{Z) « flo\Z\'^, with a < 1/2, 
Eqs. ([3T|) and ((29| yield, in the limit of /_i ^ oo: 



{p/^lo) In {^l/^lo) , for « = 1/2, 

a (1 - 2ay^ W^^o)^^" , fort a < 1/2 



(33) 

where constant fiQ is determined by the structure of 
n±{Z) at finite values of Z. Note that both results ((32|) 
and (|33p satisfy the "anti-Vakhitov-Kolokolov" (anti- 
VK) criterion, dN/d^ > 0, which implies the dynamical 
stability of the trapped modes in the case of the self- 
repulsive nonlinearity j27| . 



V. SELF-TRAPPING IN THE SELF-REPULSIVE 

BEG AND NON-INTERAGTING FERMI GAS 
DUE TO THE RADIAL MODULATION OF THE 
ID GONFINEMENT 



from the 3D Gross-Pitaevskii equation: 

= -^l^j + hLo^z)^! + 4a.|$P$, (28) 

where iL)± is the same trapping frequency as above [cf. 
Eq. (|22|)]. m the atomic mass, Os > the scatter- 
ing length of the inter-atomic interactions, and the total 
number of bosons is given by 

/ + 00 
|$(z)|^dz (29) 
-oo 



The mechanism of the formation of the localized 
ground state in the non-interacting Fermi gas subject to 
the action of the HO trapping potential, whose strength 
grows from the center, can also be applied to the case 
when the potential acts only in one direction, z, while its 
tightness, w^, is made a growing function of the trans- 
verse radial coordinate, r. Similar to the situation con- 
sidered in the previous section, this mechanism applies 
not only to fermions, but also to the self-repulsive BEG 
under the same type of the confinement. 

In terms of such a BEG setting, a system of equations 
for the 2D mean- field wave function, $ {x,y^t), and the 
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effective thickness of the condensate, r]{x, y), was derived 
inRef. [H: 



5$ 

7(r) 



1 



Vi + 7(r)77"^$ + (1/4) (f]-'^ + ri^) 



l + 7(r)|$rr,, 
2as\/2-KmuJz{r)/h, 



(35) 
(36) 



where > is the scattering length of the repulsive in- 
teractions between bosonic atoms, and is the Lapla- 
cian acting in the plane of (x, y). A straightforward anal- 
ysis demonstrates that Eqs. and ([55]) can be reduced 
to a single 2D nonpolynomial Schrodingcr equation. 



r(r)|$|4/3 



r(r) EE -(7(r))'/^ 



(37) 
(38) 



provided that the local density of the condensate is large 
enough: 

7|$P>1. (39) 

Stationary solutions to Eq. (|37p with chemical potential 
fj. are looked for as 

$(a;,y,t) = e-*^*C/(r), (40) 

with real function U{r) obeying the radial equation, 



1 f(PU 



-f^)+r(.)c/V3, 

r dr 



the respective number of atoms being 



TV = 27r 



{U{r) f rdr. 



(41) 



(42) 



Under the same conditions, a similar equation for a 
functional-density wave function U can be derived for 
the non-interacting Fermi gas by the reduction of the 3D 
density-functional description to 2D, as shown in Ref. 

M- 

MC/ = -ViC/ + C2D(c^.(r•))'/'C/'/^ (43) 

where C2D = ^3/5(6/(2s + l))^/^7r, and s is the semi- 
integer spin of the fermions. Note that, in the bosonic 
and fcrmionic settings alike, Eqs. (|37)) . ((38)) and (|43)) 
demonstrate that the coefBcient in front of the nonlinear 

1 /3 

term with total power 7/3 is proportional to (^^(r)) 

The ground state can be constructed, as above, by 
means of the TF approximation, which neglects the 
derivatives in Eq. (PT|) : 



(44) 



hence the corresponding dependence between the chem- 
ical potential and norm (j42p takes a simple form. 



TVtf = G/i3/^ G=2tt 



rdr 



[r(r)] 



3/2 ' 



(45) 



cf. Eqs. ([32]) and (|33l). Note that the TF approximation 
()44p yields a fully continuous wave function (similar to 
that obtained by means of the TF approximation in Ref. 
[isj . in the model with the spatially growing strength 
of the cubic self-defocusing nonlinearity in the nonlinear 
Schrodinger equation), unlike the ones found above in 
the form of Eqs. ([24| and ([3T|). which imply divergence of 
the derivatives at the boundary between the nonzero and 
zero parts of the solution, ujj_(z)r ~ -y/^/Iand 0^(Z) = /i, 
respectively. 

According to Eq. (|45)) . the norm of the trapped mode 
converges if r(r) grows at r — > 00 faster than r^^^, i.e., 
as it follows from Eqs. ([55]) and ([55)1 . the transverse- 
trapping frequency, uiz{r), must grow faster than r^, i.e.. 



UJz{r)/r^ 



(46) 



at r — >■ (X). Note that the similar condition obtained in 
the setting considered in the previous section was much 
weaker, requiring only that lj_i_{z) had to grow faster than 
■\/|z|at \z\ — >• 00, see Eq. ([25]) . It is also worthy to note 
that the trapped TF modes exist with ^ > 0, as well as 
solutions (|24)) and ((3T|) considered in the previous section. 

In the present notation, the condition p9p of the appli- 
cability of Eq. (|37)) for the bosonic gas takes the simple 
form, if TF approximation ()44|) is used: 



3/2 



> 1. 



(47) 



This condition guarantees the applicability of Eq. ([57)1 
not only to the core of the trapped mode but also to its 
tail, which decays at r — > 00, i.e., to the entire trapped 
mode. 

To illustrate the realization of the general setting con- 
sidered in this section, one may take the modulation func- 
tion as 



r(r)=ro {rl+r') 



2\/3/2 



(48) 



with /3 > 4/3 [this condition is necessary to meet condi- 
tion dm)]. Then, Eqs. (gl]) and (jH]) yield 



UTF{r) 



N 



Ait 



3/4 



(-0^ + 



2\ -3/3/8 



3/2 



(49) 



3/3-4 \ro< 



cf. Eq. ([5^ . Note that the latter N{ii) dependence also 
satisfies the anti-VK criterion. Finally, it is easy to check 
that, in addition to condition (|Tf)) of the applicability of 
the underlying equations ([57]) and (HI]), the condition 
for the validity of the TF approximation (|^ . i.e., the 
possibility to neglect the derivatives in Eq. (|4T|) with 
r(r) taken as per Eq. (|^5)) . amounts to /i ^ Tq"^, i.e., 
the trapping frequency lOz should not be too small at 
r = 0. 
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VI. CONCLUSIONS 

In this work, we have demonstrated that the recently 
proposed mechanism for the creation of bright sohtons in 
BEC and nonhnear optics, by means of the self-repulsive 
nonlincarity with the strength growing at r — oo, may 
be also realized in Fermi gases. For the two-component 
spin-balanced gas, this may be achieved by making the 
repulsion between the spin components accordingly mod- 
ulated in one or two directions, and the application of 
the ordinary HO (harmonic-oscillator) trapping poten- 
tial in the other direction(s). The analysis is based on the 
Euler-Lagrange equation produced by the minimization 
of the TF (Thomas-Fermi)-vW (von Wcizsacker) single- 
orbital density functional, which is quite reliable for the 
description of dilute normal Fermi gases at zero temper- 
ature [l^. We have concluded that the vW gradient 
term gives nearly negligible corrections to the TF ap- 
proximations, even for relatively small numbers of atoms. 
We have shown that both longitudinal and transverse 
widths of the localized ground state in the Fermi gas can 
be efficiently controlled by tuning the s-wave scattering 
length of repulsive interactions between spin-up and spin- 



down fcrmions, or by varying the number of atoms in the 
ground state. Further, we have demonstrated that the 
localized ground states can be created too in the non- 
interacting Fermi gas (in particular, in the spin-polarized 
one) , trapped in one of two directions by the HO poten- 
tial whose tightness grows fast enough in the remaining 
direction(s), from the center to periphery. It has been 
also demonstrated that the latter mechanism may cre- 
ate self-trapped ground states in the self-rcpulsivc BEC 
subject to the same confinement. 

The analysis reported in this paper can be extended 
further. In particular, it may be interesting to apply it 
to Bose- Fermi mixtures [l^, in the same settings which 
were studied here. 
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